Stripe structure, spectral feature and soliton gap in high T c cuprates 
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We show that the lightly doped L^-zSr^CuCU can be described in terms of a stripe magnetic 
structure or soliton picture. The internal relationship between the recent neutron observation of the 
diagonal (a;=0.05) to vertical (x > 0.06) stripe transition, which was predicted, and the concomi- 
tant metal-insulator transition is clarified by this solitonic physics. The phase diagram with the 
unidentified transition lines between antiferromagnetic to stripe phases, the doping dependence of 
the modulation period, the origin of the mid-infrared optical absorption are investigated compara- 
tively with other single layer systems: La2- a; Sr I Ni04 and (La^d^-zSr^CuCU. The novel type of 
quasi-particles and holes is fully responsible for metallic conduction and ultimately superconductiv- 
ity. 



PACS numbers:74.25.-q,75.25.+z,75.60.Ch 

Much attention has been focused on electronic proper- 
ties in underdoped high T c cuprates recently because it 
is hoped that understanding of the unusual ground state 
may lead us to a clue to a high T c mechanism. However, 
in spite of intensive theoretical and experimental works in 
the past little consensus for the normal state properties 
has emerged yet. 

Recently a remarkable series of elastic neutron exper- 
iments §-§ on Laa-zSr^CuCU (LSCO) has been per- 
formed, revealing static magnetic incommensurate (IC) 
structure: (1) For x=0.12, 0.10, 0.08 and 0.06, which are 
all metallic and superconducting, the modulation vector 
Q are characterized by Q = 5 ± e ) an d (^±£,^) in 
reciprocal lattice units (r.l.u.). The IC modulation runs 
vertically in the a (&)-axis of the CuC>2 plane. (2) For 
x=0.05 which is insulating the above four superspots ro- 
tate by 45° around the (5, ^) position, characterized by 
(i±e, i±e). The IC modulation is now in the diagonal di- 
rection. (3) The incommesurability e is given by e = x for 
0.05 < x < 0.12 beyond which e saturates (see Fig. 2). (4) 
For £=0.03 and 0.04 the elastic peaks are also observed 
at the commesurate (C) position, but the peak profiles 
are abnormally broaden, suggesting that the IC modu- 
lation is hidden. The static nature of these orderings at 
low T's are also corroborated by /iSR measurement Q|. 

These lightly doped systems were known to be in spin- 
glass region (0.02 < x < 0.05). In the further dilute 
region (0 < x < 0.02) upon lowering T the antiferromag- 
netic (AF) phase changes into a "spin-freezing" struc- 
ture at the Johnston line given by T/=815:r (K), below 
which the internal magnetic field probed by the La-NQR 
is mysteriously broaden [^J. This transition line will be 
interpreted differently below. 

Prior to these experiments, there has been known 
the static magnetic structures in superconducting (La, 
Nd^-zSr^CuCU (Nd) and insulating La2-a;Sr x Ni04 
(Ni). In the former Nd system [|| the vertical stripe with 
the same e = x relation as in LSCO for x > 0.07 and in 
the latter Ni system the diagonal stripe with e — x/2 up 
to x ~ 0.33 are observed 0]. Note that in YBCO6.6 verti- 



cal IC fluctuations are also reported in inelastic neutron 
experiments ||. 

Almost a decade ago one of the authors J|,[l(| predicts 
some of the above features of the static IC spin mod- 
ulation in lightly doped cuprate systems, stressing the 
charged stripe or the solitonic structure as a convenient 
and universal "vehicle" to accommodate excess carriers 
in an otherwise commesurate antiferromagnet within a 
mean-field treatment for a simple two-dimensional Hub- 
bard model. Some features of them were independently 
found by others |Tl[| at that time and are now confirmed 
by a more sophisticated method, such as DMRG Jl2| . 
Specifically we have predicted j^JU^]: (1) The diagonal 
to vertical transition occurs as doping proceeds. (2) 
The incommesurability e is given by e = rih/2 where 
rih is the hole density per site and proportional to x. In 
other words, once holes are introduced in Cu02 plane 
the ground state is always described by a modulated IC 
structure. (3) The newly created solitonic midgap band 
situated between the insulator main gap accommodates 
excess carriers, keeping the AF formation energy intact. 
The band width of the midgap state for the vertical stripe 
is wider than for the diagonal stripe, implying that the 
former (latter) has a tendency toward a metal (insulator). 

Encouraged by the success of our predictions, we now 
further extend it to a slightly more realistic case: namely, 
we introduce a next nearest neighbor hopping t' to de- 
scribe the metallic situation, which is certainly important 
to mimic the actual Fermi surface topology fTs]] . (The 
previous solutions with only the nearest neighbor hop- 
ping t were all insulating and fail to describe a metallic 
state |9|-|llJ|.) Thus the main purpose of this paper is to 
establish the solitonic picture to resolve various myster- 
ies or paradoxes in lightly doped regions, mainly focus- 
ing on LSCO. These include (A) why the metal-insulator 
transition and the diagonal-to-vertical transition occur 
simultaneously at x ~ 0.05, (B) why the spin freezing 
(0 < x < 0.02) below the Johnston line and the spin-glass 
regions exist, and (C) why the midgap state observed by 
the mid-infrared optical absorption J14| grows progres- 
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sively with x. Here we are going to give a self-consistent 
picture based on soliton physics, which is quite different 
from an influential idea of frustrating charge segregation 
by Emery et al |lE|| . 

We start out with the standard Hubbard model H = 
—Tii t j ta ti t jcl a Cj i(r + WSirii^n^i in two dimensions where 
tij = t for the nearest neighbor pairs (i,j) and tjj = t' 
for the next nearest neighbor pairs situated on a diagonal 
position in a square lattice. (The energy is measured by 
t and t' /t = —0.1 and U/t = 3.0 are chosen in the follow- 
ing.) The mean field {ni_ a } — rii + aMi at i-site is intro- 
duced where nj(-Mj) is the charge (spin) density. Thus 
the one-body Hamiltonian Hmf — — ^i,i,o-*ij c J a c j," + 
UT,i((ni^)rii t i + («i,x)«i,T) is solved self-consistently in 
/c-space under a given hole concentration or filling 1 — n/, 
per site (n/i=0 is just the half- filling) , assuming a spa- 
tially periodic structure for spin and charge modulations 
with e = rih/2. 

The resulting phase diagram in T vs rih is displayed in 
Fig.l where the simple commensurate AF state denoted 
by C is never stable at T=0. Upon raising T the vertical 
incommensurate (VIC) or the diagonal incommensurate 
(DIC) phase are crossed-over to the AF because of an en- 
tropy effect. These spatially modulated spin structures 
or soliton structures are characterized by a periodically 
placed 7r-shifted domain wall (or a kink) which neatly 
accommodates excess carriers. The charge carriers align 
perpendicularly to Q — (|, |) (see Fig.l in Ref. (l(J). The 
transition from DIC to VIC upon increasing can be 
understood as arising from the competitive effects be- 
tween the nesting of Q- vector (i ± e, | ± e) in DIC at 
low doping and gapping at the (1,0) saddle point posi- 
tion with large density of states in VIC. Note that the 
nesting deteriorates with rih, necessarily inducing DIC to 
VIC transition. 

The Johnston line T/ = 815.x (K) for < x < 0.02 
which indicates the anomalous broadening of the La res- 
onance line |l6| ] thus named a spin freezing T/, signals 
the change of the internal magnetic field distribution. It 
may correspond to the C-IC transition line in Fig.l. Al- 
though the present numerical calculation gives it a first 
order, actual transition is of second order (see the de- 
tailed analytical discussions in Ref. ||). According to 
Fig.l the known spin-glass region (0.02< x < 0.04) is 
also described by the DIC structure. Indeed the recent 
neutron Jl]-||] and ^SR § experiments suggest this is the 
case as mentioned above. Previous spin-glass identifica- 
tion simply reflects dirt effects. 

The incommesurability e is calculated as e = near 
half-filling at T=0. The reduced fc-space area spanned 
by this Q-vector is exactly filled by excess holes, dif- 
fering from the simple so-called 2fc^-instability in one- 
dimensional case (hp the Fermi wave number). It means 
that the charge stripe is completely filled with holes, 
never stabilizing the "half-filled" stripe in our calcula- 



tions |l7y . This is in sharp contrast with Tranquada's 
claim for Nd systems Q. The observed e is precisely 
given by e = x both for LSCO and Nd systems at low 
doping levels of x < 0.12. This implies e = rih which 
corresponds to a half-filled stripe if we identify rih = x. 
The observation e — x is also contrasted with the neutron 
result e = x/2 for Ni systems H, which does agree with 
our estimate. This "factor 2" paradox in LSCO and Nd 
is quite intriguing and may be important in evaluating 
the actual hole concentration in Cu02 plane. Here we 
suggest a possible scenario to resolve it: Fig. 2 shows the 
observed data of e for three systems (filled symbols) and 
n e f f (open ones) estimated by Uchida et al Jl4| through 
optical conductivity measurements. It is seen that ac- 
cording to Uchida's data under a given x the supplied 
n e ff in Cu02 is doubled for LSCO and Nd compared to 
Ni at least for x < 0.1. This immediately resolves the 
above paradox, namely, rih — 2x for LSCO and Nd. The 
theoretical relation e = n^/2 is valid for all three systems 
commonly and there is no half-filled stripe. It is noted, 
interesting enough, from Fig. 2 that a precise coincidence 
between e (filled circles) and n e ff (open circles) against 
x is seen to be hold even beyond x ~ 0.12 up to x ~ 0.25 
for LSCO and Nd. It means that beyond x ~ 0.12 doped 
holes do not come into the relevant magnetic band. 

The spatial magnetic structure is characterized by a 
soliton form with many higher harmonics of the Fourier 
components of the order parameters; Miq (spin) and uiq 
(charge). The period of the charge density is half that 
of the spin density. The tending limits of those toward 
rih ~ ¥ (half-filling limit) are given by M( 2 l+i)q/Mq — ► 
2^-, Mq/Mc -> f and ti 2 iq -» (I = 1, 2, 3, • • •) where 
Mc is the AF order parameter. That is, the spin mod- 
ulation becomes a squaring-up cnoidal wave form, far 
from a simple sinusoidal one and the charge modulation 
becomes irrelevant when rih — > 0. These features well 
coincide with the previous analytic solution sn(x, k) for 
t'=0 §. 

Fig. 3 shows the band structures for VIC and DIC. The 
original single band is fold back into the small reduced 
zone. In the AF case at half-filling the large main gap 
separates it into the filled and empty bands. The doping 
creates a new band or the soliton band, seen from Fig. 3, 
in between the AF gap whose magnitude is relatively in- 
tact upon doping (also see Fig. 7 in Ref. [0). In DIC the 
midgap band is empty to accommodate holes and situ- 
ated above the chemical potential fi(= 0). The midgap 
band width for VIC is substantially wider that that for 
DIC. In VIC the dispersive midgap band touches /i, mak- 
ing it metallic since the doping increases the soliton band 
width. 

In Fig. 4 the one-particle spectral weight is shown on 
the energy vs wave number plane along the symmetry 
lines for a metallic state: The reorganized band opens a 
large gap at (0,1) and {—h, \) positions where the soli- 
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tonic midgap band also situates. The hold-backed bands 
get spectral weight toward (—1,1) direction after taking a 
plateau at (0,1). This characteristic feature here is strik- 
ingly similar to that observed by photoemission experi- 
ments Q. The midgap band creates "quasi-particles" 
in its bottom at (0,1) and simultaneously the valence 
band also produces the same amount of "quasi-holes" at 
the valence band top at (— i, |). These newly created 
quasi-particles and -holes are mobile, fully responsible 
for metallic conduction and ultimately superconductivity. 
The Fermi surface pockets for these novel quasi-particles 
and holes are situated around the (0, 1) and (—5, 5) po- 
sition respectively in the Brillouin zone. The spectral 
weight has a two-fold symmetric form. The midgap band 
does not touch fi(= 0) at (1,0). 

The aligned holes in VIC on a wall move easier than 
in DIC because of the difference of the effective hoppings 
(vertical vs diagonal hopping processes), resulting in the 
different effective masses for two soliton bands. The band 
dispersion is very anisotropic; the effective mass perpen- 
dicular to the Q-vector is much heavier than that par- 
allel to the Q-vector. The holes localized on a wall are 
mobile only along the domain wall, once metallicity is 
attained upon doping. These make VIC more metallic 
while DIC insulating. In fact in our calculation in Fig.l 
VIC (DIC) just corresponds to a metal (insulator), as co- 
inciding with the observation EE]. However, we remark 
that the exact one-to-one correspondence is not generic, 
but the tendency that VIC is easier to become metallic 
than DIC is generic, independent of the choice oft'. Note 
that Ni systems are all DIC B and insulating, and Nd for 
x > 0.07 |§ and YBC0 6 . 6 g are all VIC and metallic, 
agreeing with this general rule. 

In Fig. 5 we show the optical conductivity for two dop- 
ings. As rih increases, the newly created midgap ab- 
sorption from the valence (conduction)-to-soliton band 
transition process progressively grows inside the main 
AF gap at u> ~ 1.7. Its area corresponds to the hole 
concentration. These characteristics agree with the com- 
mon optical observations for all three systems H]. The 
long-standing mystery of the origin of the mid-infrared 
absorption can be attributed to the soliton band. 

In conclusion, we have found: (1) The metal-insulator 
transition (MIT) is accompanied by the profound mag- 
netic structural change from the diagonal to vertical 
stripes as x increases. Note, however, that the actual 
MIT may be more complicated in the sense that the Q 
vector may change gradually, resulting in a second-order 
like MIT. It is interesting to carefully trace the super- 
spots in neutron experiment for the critical concentra- 
tion of ir=0.05 under varying T. (2) The incommesurate 
spin modulation persists down to x — > at T=0, even 
deeply into the known "commesurate" (0 < x < 0.02) 
and spin glass regions (0.02 < x < 0.04) under the 
assumption that the doped holes spread out the whole 
sample without any defects or impurities. We point 



out a possibility that the Johnston line may signal this 
commesurate-incommesurate transition. (3) The soli- 
tonic midgap band formed upon doping is responsible 
to the so-called mid-infrared absorption in optical con- 
ductivity measurements. (4) The incommesurability e; 
deviation from (5, 5) position is linear in the hole con- 
centration rih per site at least in lightly doped region, 
coinciding with the recent neutron observation by Endoh 
and Yamada group |l]-|i| . 

Finally we remark that judging from the optical data 
of the main insulating gap (Ni~4eV and LSC0^2eV) 
the Ni system may have larger U than in LSCO [ p^| . 
This makes DIC in Ni stabler, keeping the system insu- 
lating up to higher doping level (x > 0.33) (see Fig. 5 
in Ref. jl0| where the critical concentration of the DIC- 
VIC transition is shown to increase with x) . This may be 
one reason why heavily doped nickelates remain insulat- 
ing. As for Nd system we predict a similar 45° rotation 
of superspots in neutron experiments and concomitant 
metal-insulator transition at lower dopings. If these pre- 
dictions are indeed confirmed experimentally, the present 
soliton picture must be a proper starting point to begin 
with in order to consider a high T c mechanism. 

The authors thank many colleagues for enlightening 
discussions on the most recent updated experiments, in- 
cluding Y. Endoh, K. Yamada, S. Wakimoto, H. Kimura, 
K. Hirota, S. Uchida, A. Fujimori, A. Ino, J. Tranquada, 
M. Arai and M. Matsuda. 
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FIG. 1. Phase diagram in T vs %. C: antiferro- 
magnetic VIC:vertical, and DIC: diagonal incommensurate 
phases. VIC (DIC) is metallic (insulating) at T=0. Note 
that n h = 2x for LSCO. 
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FIG. 2. Incommensurability e vs doping x (left scale). 
e = x for LSCO §§ (filled circle) and Nd § (diamond). 
e — x/2 for Ni jjj (filled triangle). The effective hole density 
n e f f vs x (right scale) estimated by Uchida Jlij (open circle: 
LSCO, open triangle: Ni). 




FIG. 3. Dispersion relations for VIC(a) and DIC(b) along 
the perimeter of the reduced Brillouin zone for nj, = 0.02. 
The chemical potential is at zero. The midgap band is seen 
in between the AF gap. 
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FIG. 4. Contour map of the spectral weight along the 
symmetry lines in the original Brilloiun zone (nh = 0.1 and 
VIC). The original cosine-like band is hold-backed and gapped 
to produce the midgap bands and the "shadow" bands at 
around /i = 0. 
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FIG. 5. Optical conductivity (arbitrary units) vs fre- 
quency W in VIC for — 0.02 and 0.06. The midgap ab- 
sorption at around u ~ 0.8 grows as doping proceeds. 
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